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Abstract 

Generating function equation has been derived for the probability 
distribution of the number of nodes with k > outgoing lines in ran- 
domly evolving special trees defined in an earlier paper arXivxond-] 



mat/0205650| . The stochastic properties of the end-nodes (k = 0) 
have been analyzed, and it was shown that the relative variance of 
the number of end-nodes vs. time has a maximum when the evolu- 
tion is either subcritical or supercritical. On the contrary, the time 
dependence of the relative dispersion of the number of dead end-nodes 
shows a minimum at the beginning of the evolution independently of 
its type. For the sake of better understanding of the evolution dynam- 
ics the survival probability of random trees has been investigated, and 
asymptotic expressions have been derived for this probability in the 
cases of subcritical, critical and supercritical evolutions. In critical 
evolution it was shown that the probability to find the tree lifetime 
larger than x, is decreasing to zero as 1/x, if x — > oo. Approaching 
the critical state it has been found the fluctuations of the tree lifetime 
to become extremely large, and so near the critical state the average 
lifetime could be hardly used for the characterization of the process. 

PACS: 02.50.-r, 02.50.Ey, 05.40.-a 



1 Introduction 

In a previous paper |IJ we defined and analyzed random processes with con- 
tinuous time parameter describing the evolution of special trees consisting 
of living and dead nodes connected by lines. The initial state Sq of the tree 
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corresponds to a single living node called root which at the end of its life 
capable to produce new living nodes, and after that it becomes immediately 
dead. The new nodes are promptly connected to the dead node and each 
of them independently of the others can evolve further like a root. It is evi- 
dent that the random evolution of this type is nothing else than a branching 
process. 

In what follows we will use the notations applied in Therefore, the 
distribution function of the lifetime r of a living node will be denoted by 
T(t), and the probability that the number v of living nodes produced by one 
dying precursor is equal to j by fj where j 6 Z. f\ 

In order to characterize the tree evolution two non-negative integer valued 
random functions ptit) and Pdif) were introduced in |IJ. p,z{t) is the number 
of living nodes, while pd{t) is that of dead nodes at t > 0. It was mentioned 
also that the trivial equality 



must be valid with probability 1 at any t > 0. Here p e {t) is the number of 
lines in the tree at the time moment t > 0. For the generating functions 



W (f,*)=X> W (f,n)* B , es g^(t,z) = J2p (d) (t,n)z n , (1) 



DO 



?1=0 



n=0 



where 



V 




n\S } 



es p 



^(t,n)=P{p d (t) 



n\S } 



the following integral equations were derived: 





and 





where 



oo 



3=0 



Z is the set of non-negative integers. 
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If T(t) = 1 — e~^*, then the evolution becomes a Markov process. 

The nodes can be sorted into groups according to the number of outgoing 
lines. Denote by n(t, k) the number of nodes with k > outgoing lines at 
the time instant t > 0. A node not having outgoing line is called end-node. 
It is obvious that an end-node could be either live or dead. Therefore, the 
number of end-nodes fi(t, 0) can be written as a sum of numbers of living 
and dead end-nodes, i.e. 

fi(t,0)=fi t (t,0)+fi i (t,0). 

Since all living nodes are end-nodes fjbi(t, 0) can be replaced by (ii(t). The 
total number of dead nodes /x<f(i) is given by 

oo 

Vd(t) = fi>d{t, k). 

k=0 

In what follows we will calculate the probability distribution of Hd{t, k) 
and investigate the properties of end-nodes 0) playing important role in 
random tree evolution. 

In order to have a deeper insight into the dynamics of the evolution 
process, we will derive an important equation determining the probability 
distribution function of the tree lifetime. 



2 Generating functions 



2.1 Distribution of v 

The basic properties of the probability distribution of the number of nodes 
in a random tree are depending mainly on the distribution law of the number 
v of living nodes produced by one dying precursor. In the sequel we will use 
the notations 



EM 

introduced already in 



qi and D 2 {z/} = q 2 + qi - q\ 
| for the expectation and the variance of v where 
'd?q(z) 



dzi 



J 



1,2, 



2 = 1 
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are the factorial moments of v. It was shown in |1J that the time dependence 
of the random evolution is determined almost completely by the expectation 
value q\. The evolution is called subcritical if q\ < 1, critical if q\ — 1 and 
supercritical if q\ > 1. In the further considerations we are going to use four 
simple distributions for the random variable v. 



2.1.1 Arbitrary distribution 

It has been shown in ]]J that the equations derived for the first and the second 
moments of fii(t) , fi d (t) , . . . are true for any distribution of v provided that 
the moments q\ and q2 are finite. This type of distributions of v is called 
arbitrary and will be denoted by a. 

Many times it is expedient to assume distributions of v to be completely 
determined by one or two parameters. 

2.1.2 Geometric and Poisson distributions 

As known the geometric and Poisson distributions are containing one param- 
eter only, and so we have 



V{u = j} 



1 I _qi 

1+91 \ 1+91 



91 1 if v e P , 



where g and p refer to the geometric and the Poisson distributions, respec- 
tively. For the sake of completeness we write 



i 



if v e g, 



q(z) 



l+(l-z)gi' 
e -(l- Z )gi ; i f z/G p ! 



and 



qi, if v e g, ( <?i(i + gi), if v eg, 

E{z/} = { while D 2 {z/} = <^ 

q x , if v G p, I q u if v G p. 
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2.1.3 Truncated arbitrary distribution 

In this case the possible values of the random variable v are 0, 1 and 2 with 
probabilities / , f\ and / 2 , respectively. This distribution will be denoted by 
t. The corresponding generating function of v is given by 

q(z) = fo + fiz + fcz 2 = 1 + qi(z - 1) + ±q 2 (z - l) 2 . 

This choice of q(z) has a great advantage, it makes possible to obtain exact 
solutions of the generating function equations of fJ>e(t) , ^d{t) , ■ ■ ■ characteriz- 
ing the tree evolution. 




0.5 1 1.5 2 2.5 

qi 

Figure 1: The permitted values of qi and qi in the case of distribution t of v. 
It seems to be useful to cite the following trivial relations: 

fo = 1 - qi + ^ q 2 , fi = qi- <?2, h = \ q2 > 
which follow from the equations 

/o + /i + / 2 = l, fi + 2f 2 = qi, 2f 2 = q 2 , 



5 



where q 2 = D 2 {^} — <?i + qf- Since f , fx and / 2 are non-negative, smaller 
than 1, real numbers and their sum is equal to 1, the possible values of q± 
and q2 are restricted. The permitted values of q\ and q 2 are shown in Fig. [I]. 
(See the shaded triangle!) 

2.2 Distribution of the number of dead nodes with k 
outgoing lines 

As has been already mentioned, a living node may create k > new living 
nodes and after that it becomes immediately dead. This node is called dead 
node of out-degree k. It was introduced the random function fid(t, k) giving 
the number of these dead nodes at time instant t > 0. Now, we want to 
determine the probability that [J>d(t, k) is equal to n > provided that at 
t = the tree was in the state Sq. First of all, we define the probability 
generating function 

oo 

gi d \t,z) = J2p ( *\t,n)z n , (4) 

n=0 

where 

p[ d \t } n)=V{fi d (t,k)=n\S } 

is the probability that at time t > the tree has exactly n dead nodes with 
k > outgoing lines provided that at t — it was in its initial state So. By 
using similar considerations as we did in 0, we have 

pf(t,n)=e- Qt 5 nfi + 



/oo ( °° 

e _ Q(t _ t0 | /q [6(jk6ni + (1 _ 6ok)6no] +^/ i ^(t',n)| dt', 



(5) 



where 



niH h«j=«— 1 i=l niH hnj=n i=l 
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The expression in square brackets at / reflects that two mutually excluding 
possibilities exist depending on whether k = or k ^ 0. One can see 
immediately that the generating function defined by Eq. (^) satisfies the 
equation 

g { k d \t, z) = e~ Qt + Q f e-^'Vofl - (1 - z) 4,o] dt'+ 

Jo 

oo 



q / e -Q(«') 
By rearranging the right side we find 



*) = e"«* - (1 - z)f 5 kfi Q T dt>- 

Jo 



\ f k (1 - z) (1 - 5 0>k ) gf\t',z) "-q gf\t\z) \ dt> . (6) 



Q 



When k > then the equation (|6]) takes the form: 



z = e 



-Qt_ 



Q / f 

'0 



The differential equation equivalent to (El) is nothing else than 



'^(f,*)! dt>. (7) 



Qg^(t,z)-Qf (l-z)6 kt0 - 



dt 



f k (1 - z) (1 - 5 , k ) g { k d \t,z) -q g k a >(t,z) 



and the initial condition is given by lim^ gj® (t, z) = 1. 



(8) 
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2.3 Joint distribution of the numbers of two dead 
nodes with different out-degrees 

The following step in our considerations is the determination of the joint 
probability distribution 



V{p d (t,ki) =n l7 fi d (t,k 2 ) = n 2 \S } = p k d ^ k2 (t,n l7 n 2 ) 



(9) 



where k\ ^ k 2 . It is clear from the definition (|J) that P k ^ k2 (t, Hi, n 2 ) is the 
probability that in the time interval (0,t) the evolution produces n\ nodes 
with k± and n 2 nodes with k 2 outgoing lines provided that at the moment 
t — the tree was in its initial state Sq. By using similar arguments applied 
in deriving the backward equation (ffi), we obtain for the generating function 



oo oo 



g^ k2 (t, Zl ,z 2 ) = E&^^l *T z ? 

n±=0 U2=0 

the following equation when ki ^ k 2 : 



(10) 



g%> M (t, 21, z 2 ) = e- Qt -Q f f e - Q ^[(l - z x ) 8 kl , + (1 - z 2 ) 5 k2 , ] dt'- 

Jo 



Q 
Q 



dt'- 



-^\l-z 2 )f k2 gZ 2 (t',z u z 2 ) 



lk 2 



dt'+ 



Q 



-Qtt-t'lnln^ (+' 



q [g k a > k2 (t', Zl , Z2 )}dt'. 



(11) 



It follows from this equation that 

M) u „ „a au(l 



Mi 



For the sake of the completeness we are giving here the differential equation 
equivalent to (|TT|). It has the form: 



Q 9 { k % 2 (t, zi, z 2 ) + Q q{g^ M (t, z u z 2 )}- 



(d) 



dt 
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Q /o[(l - Zl) <W + (1 - 2&) <J*a,o]- 



(1 - ^l) /fci [#2jfc 2 0, 21, 3j)J - (1 - Z 2 ) f k2 g^ M (t, Zi, Z 2 ) 

with the initial condition 

^^ fc2 (t, 21,22] = 1- 



fei 



k 2 



(12) 



2.4 Distribution of the number of end-nodes 

The dynamics of the random tree evolution is controlled by the the end-nodes. 
It is evident that the living end-nodes are responsible for the development of 
a tree, while the dead end-nodes represent those points where the evolution 
was stopped. The probability distribution of the number of dead end-nodes 

Ad) 



p (t, n) can be obtained by substitution k = into Eq. 
show that the generating function 



It is easy to 



9 { d \t,z) = Y,Po d \t,n) z" 



n=0 



satisfies the equation 



^,z) = e-«*-(l-z)/ (l-e-« t ) + Q / e 



-Q(t-t') 



Q 



gf\t\z)} dt'. (13) 



In order to have an insight into the interplay between the living and dead 
end-nodes it seems to be useful to calculate the probability distribution of 
the random function 

/i(t, 0) = fi t (t, 0) + fi d (t, 0) 



and the joint distribution of fJ>i(t, 0) and [id(t, 0). 

With help of similar arguments used for the derivation of Eq. 
be easily obtained the equation 



it can 



g (t,z) = e' Qt z -(l-z)f (l-e 



-Qt- 



+ Q I e-^'hW.z)} dt', (14) 
'0 
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for the generating function 

OO CO 

g (t,z) = 0) = n\S } z n = £p (*,ra) z n , 



n=0 



n=0 



where po(t, n) is the probability that the number of all end-nodes at t > is 
equal to n provided that at t — the tree was in the state <So- 

Now we would like to determine the joint distribution of fj,e(t, 0) and 
fid(t,0), i.e. the probability 

V{ne(t) =n 1 , fi d {t,0) = n 2 \S } = Po'^it, n u n 2 ). 
It is evident that 



py \t,n 1 ,n 2 ) =e Qt S niA 5 n2fi + 



+Q 



f 

Jo 



-Q(t-f) 



/o<5 ni ,o<W + £/fc R e ' d \t',n 1 ,n 2 ) 



k=i 



dt', 



where 



R^ d \t',n l ,n 2 ) = E U^V,"^). 

rinH hn lfc =ni 7i2lH hri2fc=™2 j=l 

Simple calculations show that the generating function 



OO OO 



g^ d \t, Zl ,z 2 ) = J2P° ,d) ^ n ^ z " * 

ni=0 ri2=0 



112 
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satisfies the equation 

gtf' d \t,z 1 ,z 2 ) = e~ Qt Zl - / (1 - z 2 )(l - 



+Q 



e -Q(t-t>) q 



(t,zi,z 2 ) dt', 



(15) 



which will be used in the next section for the determination of the correlation 
between the random variables ^{t, 0) and fi d (t,0). 

Finally, we would like to cite the result of calculations concerning the 
probability to find the random function 



k=i 
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i.e. the number of non end-nodes at t > to be equal to n, provided that at 
t = the tree was in the state So- Denote this probability by p d °\t,n) and 
let us introduce the generating function 

oo oo 



n=0 



n=0 



It can be proven that it satisfies the integral equation 



gf{t,z) = e-^+f (l-z)(l-e-^)+zQ / e 



-Q(t-t') 



g d °\t',z)} dt', (16) 



which is equivalent to the differential equation 

dgf{t,z) 



Qfo(l-z)-g^(t,z)+Qzq g ( d °>(t,z) 



dt 



with initial condition \im t _, g^\t, z) = 1. 



2.5 Average characteristics 

It is difficult to find solutions of the generating function equations even in 
those cases when the distribution of v is known and simple. Therefore, 
in this section we would like to deal with the average properties of tree 
evolution, and will derive equations for the expectation values and variances 
of the number of nodes of different kind. In order to have an insight into 
the character of the stochastic interplay between the numbers of the living 
and dead end-nodes, we will investigate the time variation of the correlation 
between these nodes. 



2.5.1 Dead nodes with k outgoing lines 

Let the first task be the calculation of the time dependence of the expectation 
value and variance of dead nodes with k > outgoing lines. By using the 
relation 

= m[ d \t,k), (17) 

2=1 

2 As has been already mentioned the truncated (t) arbitrary distribution of v is one of 
the rare exceptions when the generating function equation can be exactly solved. 



E{fj, d (t, k)} 



d_g%z) 
dz 
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from the probability generating function (pi) we obtain that 



m [ d \t,k) = f k (l-e' Qt ) + qi Q / e- Q ^ m[ d \t',k)dt>. (18) 



The solution of this equation can be written in the form: 

[1 - e -( 1 -9i)Q*] , if qi ^ 1, 



m[ d \t, k) 



1-91 
fk Qt, 



(19) 



if qi = 1. 



Since 



^/i d (t, fc) = fi d (t), 



k=0 



it is obvious that 

oo 

J2 m i\t,k) = mf\t) 



k=0 



_J_ [i_ e -d-«)0t] j if 9l ^l, 
Qt, if g! = 1, 



what has been already derived in |]. By using this relation we can conclude 
that 



m[ d) (t, k) 
mf\t) 



fk 



Vk e z. 



(20) 



In order to calculate the variance D 2 {/^(£, k)} we need the second facto- 
rial moment m^(t, k). It can be shown that m£\t, k) satisfies the integral 
equation 

m f(t, k) = q x Q [ e-W"^ mf (*', k) dt'+ 



Q 



the solution of which can be written in the form: 



\2kf k m[ d \t',k) + q 2 \m (d \t',k) 



dt' 



m (d \t,k) 



fk 



1-91 



2k + 



Q2 



1-qi 



:i-e- at )- 



fl 



l-qi 



2k + 



92 



l-9i 



Qt e~ at + 



ff 



;i-9i 



q 2 i (1 - < 



■at\ 



(21) 
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for all gi ^ 1. Here the notation 

a = {l-q 1 )Q 
has been used. If q\ = 1, then the solution has the form: 

m {d \t,k) = f 2 (Qt) 2 (k + lq 2 Qt). 



Taking into account that 

D 2 {/i d (t, k)} = m { 2 d \t, k) + m[ d \t, k) [l - m[ d \t, fc)] , 

where 



fk 



1 - 



m^A:) l-m^A;) 



i-9i 



(l-e" Qi ) + 



e~ at (l-e~ at ), 



we obtain the following expression for the variance: 

B 2 {fi d (t,k)} = 



fk 



1-91 



1 + 



fk 



l-9i 



l-9i 

fk 



2k - 1 + 
2A- 



92 



l-9i 

92 \ /4J -„/ 



(l-e- Qt ) 



(l-9i) ; 



1 + 



l-9i 



e- at (l -e~ ai ), 



l-9i 

92 \ - n l i -| 



(22) 



(23) 



9i 7^ 1 and VA; e 2, 

where g 2 can be replaced by D 2 {z/} — q 1 (1 — qi). When the evolution is 
critical, i.e. when q 1 — 1, then we have 

V 2 {^ d (t,k)} = f k Qt 



l + (k-l)f k Qt + -f k B 2 {u} (Qt) 2 



(24) 
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If the time t is converging to oo the Eqs. ( P3|) and ( |2~4]) show that the 
variance remains finite in the case of subcritical evolution only. Introducing 
the notation 

q 2 D» + (1 - 9l ) 2 



we obtain immediately the formula 



l-qi 



fk 



lim T> 2 {fi d (t,k)} 

t— >oo 



I — ( 2(k — 1) + <b 



1 - qi \ 1 - ?i 

if ?! < 1. 



Vfc G Z, 



(25) 



2.5.2 Properties of end-nodes 

Now we would like to deal with some properties of end-nodes. The time 
dependence of the expectation of the number of end-nodes is one of the 
simplest characteristics of randomly evolving trees. It can be calculated 



from Eq. (|14]). By using the relation 

E{p(*,0)} = mi(*,0) 



dg (t,z) 
dz 



2 = 1 



after some elementary mathematics we have 

Jo_ + hJpqi e -(l-?l)Qt ; if qi ^ 1; 



mi(t,0) 




(26) 



if 9l = l. 



It is interesting to note that the expectation value of the number of dead 
end-nodes m[ (t, 0) can be easily calculated from (p6|). By substituting 



ni 



f\t,0) = e - ( - 1 ~ qi)Qt into the equation 



mS d) (t, 0) = mi (t, 0) - r< j (t , 0) , 



one gets the formula 
mf°(*,0) 




,-(l- qi )Qt 



), if 9i ^ 1, 



if 9i = 1 



(27) 
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1 2 3 4 5 

time in Q" 1 units 



Figure 2: Time dependence of the ratio of the expected number of end-nodes mi(i,0) 
to that of all the nodes mi(t) in subcritical (q\ — 0.9), critical (q% = 1) and supercritical 
(qi = 1.1) evolutions. 



which is the same what follows from Eq. ([19]) when k = 0. 

In many cases it seems to be important to know how the evolution process 
does alter the ratio of the expected number of end-nodes mi(t, 0) to that of 
all the nodes mi(t). Taking into account that 



mi(t) 



1 + Qt, if qi = 1, 

after elementary calculations we have 

/o+(i-/o-gi) e-d-giw -r /, 
mi(t,0) J e-d-n)«* ' 11 



mi(t) 



i+/o Qt 



if gi = 1, 



5 See Eq. (45) in 
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and Fig. |2] shows the ratio vs. time plots for values q\ = 0.9, 1, 1.1 at 
/o = 0.51. As seen the curves are approaching the limit values 



mi (t,0) 
hm -— 

t^oo mi(t) 



fo, 



if 9i < 1, 



1 - if gi > 1 



quite rapidly.Q 

In order to calculate the variance of the number of end-nodes we need 
the second factorial moment of fj,(t, 0). By using the relation 



v d 2 g (t,z) 



dz 2 



J 2 = 1 



m 2 (t,0) = E{//(t,0) ^(t,0)-l]} = 

we obtain from Eq. ([14]) the integral equation 

m 2 (t,0) = gi Q [ t e- Q ^m 2 (t',0)dt' + q 2 Q P e ~QW [ mi (t',0)} 2 dt 1 , 
Jo Jo 

the solution of that can be written in the form 

2 -i —at 



m 2 (t,0) = q 2 



o 



1 - e 



1 - 9i J 1 - 9i 



+ 2q 2 /o( )~ /o 7 2 9l) Qt e-«*+ 
(1 - 9i) 2 



+92 



1 - /o - ?1 



-a* 



1-e 



-at 



l-9i 



l-9i 

if 9l ^ 1. 

When qi = 1, i.e. when the evolution is critical we get 



Qt + fo (Qt) 2 + 3/0 (^) 3 



m 2 (t,0) = q 2 

Finally we have the variance of //(i, 0) in the following form: 



D 2 {M*,0)} 



foil ~fo- Qi] 
(l-9i) 2 



92 



1-qi V 1 - 9i 



fo 



'l-e- at ) + 



4 It is to note that 1 - > f if q 1 > 1. 
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2 q2 fo{ )- f °- qi) Qt e-*+ 

(i -qiY 

2 



1 + f 1 ~ /o ~ gl ] e- Qt (l-e-^), (28) 



l - 91 / V 1 - 9i 

if 9i ^ 1, 

and 

D 2 {/i(t, 0)} = (q 2 - /„) (1 + / Qt) Qi + ^g 2 / 2 (Qt) 3 , (29) 
if gi = 1. 

It can be easily proven that if q\ = 1, then (72 — /o > 0- Since gi = X^fcli = 1 an( i 
J2kLo f k = 1 it is evident that 



and hence 



. k, 

k=2 k=2 



/o = X>-l)/ fc . 

k=2 

By using this expression of /o we find that 

00 00 
<?2 - /o = $>(fc -!)-(*- 1)] A = J> - I) 2 A > 0. Q.E.D. 

k=2 k=2 



As follows from Eqs. (|28|) and ( f2"9|) when t — > 00 then the variance of 



fi(t, 0) tends to finite limit value in the case of subcritical evolution only, 

In order to show the main features of the time variation of fluctuations 
occurring in the end-node number of trees the relative variance and the rel- 
ative dispersion of /i(t,0) vs. time have been calculated. The results are 
plotted in Fig. £| It has been assumed that the distribution of v is arbitrary 
in that sense as it was defined in section 2.1.1. In the upper part of the Fig. | 
one can see that the relative variance of fj,(t, 0) reaches a maximum just after 
the beginning of the process but in the case of subcritical and supercritical 
evolutions only. If t — *■ 00, then we have 



lim 



D 2 {/i(t,0)} fo + 11 qi < X ' 



ooE 2 {/iM)} 



92 

gi-i 



if 9i > 1, 
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Figure 3: Relative variance (upper) and dispersion (lower) of the number of end-nodes 
vs. time in subcritical, critical and supercritical trees. The distribution of v is arbitrary 
and D 2 M = 0.9. 

The proof of the inequality q 2 > q\ — 1 is very simple. If we substitute qi and q\ by 
Y^k=\ — 1) fk an d SfcLi k /fci respectively, then we can write 

oo oo 

£ fc(fc - 1) / fc - J2 k fk + 1 > o, 

k=l k=l 

i.e. 

oo 

X> - i) 2 A + /o > o, 

and this is trivially true. 
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Figure 4: Relative dispersion of the number of end-nodes in subcritical trees vs. time. 
Curves a, g, and p correspond to the arbitrary, geometric and Poisson distributions of u, 
respectively. In the case of arbitrary distribution: D 2 {^} = 0.9. Upper: < Qt < 100. 
Lower: < Qt < 5. 



If the evolution is critical 
D 2 {M*,0)} 



1), then the relative variance 



EHM^O)} 



(?2 - fo) 



Qt 



1 + foQt 3 



(Qty 



i + foQt) 2 



is increasing monotonously to infinity with t. The curves in the lower part 
of Fig. |3| show the time variation of the relative dispersion of /i(t,0) which 
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has a finite limit value 



t-~> E{ M ((,0)} 1-9, \1 -9i 

in the the subcritical evolution only. 

The influence of the distribution law of 1/ on the relative dispersion can 
be seen in Fig. [|. The calculations have been carried out in the case of 
subcritical evolution. In the upper part of the figure it can be seen that each 
curves tends to a finite asymptotic value when t =>- oo. It is remarkable 
that the geometric distribution of v brings about much larger fluctuations 
in /i(t, 0) than the arbitrary and Poisson distributions. The lower part of 
the figure shows the beginning of the time dependence which is reflecting the 
effect of two competing processes. One of them is the formation of new living 
nodes, while the other one is the death of end-nodes. 

It seems to be worthwhile to calculate the time dependence of the expec- 
tation and variance of the number of dead end-nodes 0). Substituting 
k = into the Eqs. (fL9) and (|23l) it can be derived the relative variance 



D 2 {/i d (t,0)} ( q 2 | 1-q, t \ 1 



E 2 {fi d (t, 0)} fo 



and 



1 - e~ at ) 2 \ 1 - q x J 1 - er at ! 
if 9i 1, 



0)}_1 + 1 



E 2 {^(t,0)} 3" z " 

if gi = 1. 

The relative variance of Hd{t, 0) converges to finite value, when t — > oo, 
in both subcritical and supercritical evolutions but diverges if the evolution 
is critical. Let us introduce the notations 



d 2 {M*,Q)} 

t-oo E 2 {/i d (t,0)} 



rt>^ , it z/ G a, 
nA , if v G g, 



ri/f\ if f G p, 
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and summarize the limit values of relative variances. We find that 

— - 1 - gi + -j^, it gi < 1, 

if gi = 1, 

if q 1 > 1. 



(a) 



OO, 



31-1 ' 



If the distribution of v is geometric, then 



TV 



(<?) 



gf i±2i. 

*1 l-qi 



OO, 



if gi < 1, 
if qx = 1, 



-^L _ i if qi > i, 

V 31-1 ' ^ 1 ' 

while if it is of Poisson type, then 

' (1-91)^ + %^, ifgi<l, 

if gi = 1, 



rv 



(p) 



OO, 



31-1 x > 



if gi > 1. 



Finally, it seems to be interesting to look at the time dependence of the 
relative dispersion the number of dead end-nodes. It follows from Eqs. (|19D 
and (EST) that 



D 2 {M*,0)} /q 
E{^(t,0)} 



92 



/o 



1 - gi V 1 - gi 



92 



+ 1 e 



1-gi \l-9i 

92 



1 + 



-Of 



2/o 



Qt 



if 9i ^ 1, 



and 



D 2 {Mt,o)} 

E{// d (t,0)} 



if 9i = 1. 



1 - gi 1 - e 



1 



-at 
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Figure 5: Relative dispersion of the number of dead end-nodes vs. time in subcritical 
trees. Curves a, g, and p correspond to the arbitrary, geometric and Poisson distributions 
of v, respectively. In the case of arbitrary distribution: D 2 jV} = 0.9 and /o = 0.51. 
Upper: < Qt < 100. Lower: < Qt < 5. 

Fig. [| shows the time dependence of the relative dispersion of the dead end- 
nodes fj,d(t, 0) in subcritical evolution for all the three (a, g, p) distributions 
of v. In the lower part of the figure one can see that the relative dispersion 
curves have minimum just after the beginning of the process. (In the case of 
the arbitrary distribution of v the minimum sites are given in Table I.) 

Table I: Sites of the minimum in relative dispersion vs. time curves at several q\ and 
D 2 {V} values. (Time is measured in Q^ 1 units and v G a.) 
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d 2 MW 


0.90 


0.95 


1.00 


1.05 


1.10 


0.90 


1.7724 


1.7223 


1.6617 


1.6069 


1.5441 


1.00 


1.5849 


1.5448 


1.5000 


1.4514 


1.3999 


1.10 


1.4332 


1.4006 


1.3636 


1.3233 


1.2803 



3 Lifetime of trees 
3.1 General considerations 

It is obvious that the evolution of a random tree will stop at that time instant 
9 which satisfies with probability one the equation pe(9) = 0. The random 
variable 9 is called lifetime of the tree. In order to determine its distribution 
function 

V{9 < t\S Q } = L(t), (30) 
one has to recognize that the probability 

V{fi £ (t) = 0\S }=p {E) (t,0) 

to find zero living node at time moment t > in a tree is the same as the 
probability that the lifetime 9 of that tree is not larger than t > 0, therefore, 
one can write 



V{6<t\S } = V{nt(t) = 0\So}, 

i.e., 

L(t) =p®(t,0) =lim g {i \t,z). (31) 

It is clear that if < t\ < t% then L(t\) < Lfo), i.e., L(t) is a non- decreasing 
function of its argument, hence the limit relation 

max L(t) = lim L(t) = < 1 (32) 

0<t<oo t^oo 

must be valid. We will call the quantity dying- out-probability, and prove 
the following theorem: 

Theorem 1 If q± < 1, i.e., the random evolution is not supercritical, then 
= I? while if q\ > 1, i.e., the evolution is supercritical, then is 
equal to the single, smaller than 1, non-negative root of the function ip{y) = 
g(y)-y, ye[0,l].B 

5 It is evident that ^(1) = 0. 
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Proof. For the proof we exploit the fundamental property of the gener- 
ating function g^(t, z) which is expressed by the equation 

g«>(t + u,z)=g«>[t,g«>(u,z)]. 



Applying the relation (31) we have 

L(t + u)=g®[t,L( U )] 

and since 

lim L(u) = lim L(t + u) = L 

u^oo u— >oo 

we can write for every t > that 

L 00 =g®(t,L 00 ). 

Putting g^(t, Loo) into 

dg«>(t,z) 
alt 

derived from (|2[), we obtain 



Qq[g®{t,z)]-Qg®(t,z) (33) 



q (L^) - L^ = 0. (34) 

Considering that q(y) is probability generating function, i.e., lim^i q(y) = 1, 
according to a well-known theorem of generating functions , it is clear that 
if q\ > 1 then Eq. (|34]) — besides the trivial fixed-point 1 — must have an 
other, smaller than 1, non-negative fixed-point Loo too, and this is what we 
wanted to prove. Q.E.D. 

Let us introduce the probability 

S(x)=l-L(x), (35) 

to find the tree at the time instant x = Qt in living state. The S(x) will be 
called survival probability. By taking into account the properties of L(x) one 
obtains 

f 0, ifgi<l, 
Q lim^(x) = <^ (36) 
I Soo = 1 - Loo, if q\ > 1. 



3 The proof of this and the following theorems is based on Sjewastjanow's ideas 
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By using Eqs. (^), ([H]), and (|3"5] ) one gets 

dS 



(h -q(l-S)-S + l, (37) 



which has the solution 



X ~l S{x) <1{1-V)+V-1' (38) 

if the initial condition is S(0) = 1. Now, we would like to derive asymptotic 
expressions of S(x) for large x = Qt in the cases of subcritical, critical and 
supercritical evolution. 

Theorem 2 If the integral 

f 1 g(i - y) + qiv - 1 , , ~ , qm 

/ ~T7i \~i TT ^ = ~ lo S K 39 

finite, then in the case of subcritical evolution the survival probability S(x) 
has the following form: 

S{x) = K e - {1 ~ qi)x [1 + o(l)] (40) 

when x = Qt =>- oo. 

Proof. The proof is simple: the identity 

e -{l-qi)x 

log = (gi - 1) ar- log5(x), 

by using the Eq. fl38|), can be rewritten in the form 

i e ~ [1 ~ 9l)x < u f 1 d y f 1 d y 

log q/ v = (gi - 1) / -7: r— 7+/ — = 

1 q(y-l)+ qi y-l 

s{x) yW-y)+y-i] y ~ { h 

and hence 

S(x) = e~ k{x) e - {1 ~ qi)x . 
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From this we find 

S{x) 



Jim = e~ k ^ 

x— »oo e 



(l-gi)a; ' 



where 

Jo vm 1 - v) + y - i] 

since S(oo) — if q± < 1. According to the assumption (|39| ) we can write 

S(x) = e- fc(oo) e - {l ~ qi)x [1 + o(l)] = if e- (1 - gi)1 ' [1 + o(l)], 

and this is what we wanted to prove. Q.E.D. 

It is to mention that in the case of t distribution of v from Eq. (^) we 
obtain 

K 



1+1 _22_' 



and so 

e -(l-9i)a; 



1 + 



2 



Theorem 3 If q 2 < oo and gi = 1, t/ien £/ie asymptotic expression for the 
survival probability is given by 

S(x) = — [1 + (1)] (41) 

when x = Qt =>- oo. 

Proof. By using the series expansion theorem according to which 

q(l-S) = l-S + ±q"[b(x)) S 2 , 



where 1 — S < b{x) < 1, from Eq. (|3T| ) we obtain 

Since if x ^ oo, then S(x) =>- we can write 

g"[6(x)] = g 2 + e(x), where lim e(x) = 0. 
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Taking into account this expression for q"[b(x)] we have 
the solution of which can be written in the form 



S(x) 



1 1 
- q 2 x H — 

2 H 2 



e(u) alv + C 



The initial condition 5(0) = 1 results in C = 1, and therefore 

-l 



S(aO = 



g 2 a; 



2 1 f 

H 1 / e(u) cfe 

g 2 x Jo 



By applying the L 'Hospital rule we find 



lim — 

x— +oo a; 



e(f ) rff = lim e(x) = 0, 

x-^oo 



and hence 



g 2 x 



what we wanted to prove. Q.E.D. 

The third task is to derive the asymptotic expression for S(x) in the case 
of supercritical evolution. 

Theorem 4 If the integral 



1 



g(i - y) + y - 1 - (gi - 1) (y - 
(j/-5oo)[g(i-y) + j/-i] 



= r(oo) 



(42) 



is /mzte in £/ie case of q 1 > 1, then 

S(x) = Sec + (1 - Soo) e 



-r(oo) -(qi-l)rr 



[1 + 0(1)], 



(43) 



/or x = =>• oo where = 1 — < 1 is the limit value of the survival 
probability. 
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Proof. From the identity 



e -(<2i-l)z 

log ^j— — = -(qi - 1) x - \og[S(x) - S 



we obtain 



e -( qi -i)x f dy f dy 

•^sw^r-^- 1 ) J 9( i-»)+v -i + J t 

8(x) S(x)-S 00 



and since 



dy [ dy 1 
+ log 



y J y-Soo s(x)-s, 



oo 



we can rewrite Eq. (03) in the form 



1 — <? 

W — °° e -(n-i)x = r ( x \ 

where 

i 

, s _ f q(l-y)+y-l-(gi-l) (y-gpo) 

S(x) 

and so we have 

S(x) =S 00 + (1- 5oo) e~ T(x) e~^~ l)x . 
It has been assumed r(oo) to be finite and so 

g(g) - gpo _ „-r(oo) 

J™ e -(«-i)« bca)e <QO - 

Thus the theorem is proved. Q.E.D. 

In the case of t distribution of v we find immediately that 

Soo = 2— -, and r(oo) = -logS^, 

and hence obtain 

S(*) = 2^ + 2^ (l-2^-) e-fe- 1 )' [1 + 0(1)], 

when x = Qt => co. 
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3.2 Some exactly solvable models 



In order to demonstrate the characteristic features of the lifetime of random 
trees we will use such probabilities fk, k = 0, 1, . . . for the offspring num- 
bers v that make possible to solve exactly the equation fl3"T|) of the survival 
probability. In the following, two special cases will be investigated, namely 



q(z) 



1 + qi(z - 1) + \ q 2 (z - l) 2 , model t, 



[1 + 9i(l " z)]- 1 



(45) 



model g. 



It is clear that the first case corresponds to the zero-one-two, while the second 
to the geometric distribution of the offspring number v. 



3.2.1 Survival probabilities 

By using Eq. fl3"7|) in the model t we obtain 

f = -il-q,)S-\ q2 S\ 



(46) 



and taking into account the initial condition S(0) = 1 we have the solution 
in the form 



S(x) 



-{l-q 1 )x 



1 + 2(T^) (! - e ~ (1 ~" l)x ) > if ?i<! 



2+<j 2 x ■ 



2 21^1 

<?2 



1 + (1 _ 2 21=1) e -(qx-l)l 



ill 



-1 



if gi = 1, 



if qi > 1. 



(47) 



The corresponding equation in the model g (i.e. when the distribution of 
v is geometric) is given by 



dS_ 

dx 



S + l = -S 



l-gi(l-S) 



1 + qiS 1 + qiS ' 

the solution of which can be easily obtained in inverse form 

ji- logi + j^r log[l- ft (l-5)], ifgi^l, 



(48) 



^ + logi 



(49) 



if gi = 1, 
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satisfying the initial condition S(0) = 1. 

The density function of the lifetime measured in Q _1 units can be calcu- 
lated by using the relation 



£(x) 



dL(x) dS{x) 



dx 



dx 



It is elementary to show that the density function is decreasing monotonously 
from £(0) = /q to zero. In Fig. 6 one can see the density function curves 
versus time x = Qt for subcritical, critical and supercritical trees and for 
both distributions of v. 

Let us calculate now the characteristic function of the random variable 
d = Q8, i.e. of the tree lifetime measured in Q _1 units. Since the moments 
j = 1, 2, ... do not exist if the qi > 1, the calculations are restricted 
to the case when qi < 1. One can write 



<p(u) = E{e-^} 



dL(x) 



x dS(x) 



-ujx(y) 



dy, 



where u is a complex number with §lu) > 0. Performing the substitution 

1 



x(y) 



logy 
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y + 2 



in the model t one obtains the characteristic function 



(p t (u) = (1 + 7) 



ui/3 



f 


y 


/o 


.2/ + 7. 



w/3 



(50) 



where 



= (i-(&; 



and 



7 



?2 



gi < 1. 



In the model g we should substitute for x the expression given by 



i.e. 



l-qi y l-qi 



log[l - gi(l - ?/)], 



and we find 



= / (1 - y) ufi (1 - fty)"^ dy, Ql < 1. 
Jo 



(51) 
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Figure 6: Time dependence of the density function £(x) of the tree lifetime $ in models 
t and g. 

3.2.2 Expectation and variance of the lifetime 

From the characteristic functions y>t(u;) and <p 9 {uj) it can be easily calculated 
both the expectation value and the variance of the lifetime $ of a subcritical 
tree. In the model t for the expectation value one has 

V du y w=0 j 1/(1 + 7) 

/? 7 log(l + -) =- logfl + ^-^-V (52) 
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while in the model g one finds expression 

W = =1-108(1-30. (53) 




Figure 7: 
and g. 



Dependence of the average value of the tree lifetime $ on qi < 1 in models t 



The curves in Fig. ^ show the dependence of the expectation of the tree 
lifetime on the average branching parameter q± in both models t and g. One 
can observe the difference between the two curves to be unimportant. In both 
cases if q\ =>- 1 then the expectation value becomes infinite like log(l — q%)~ x - 

For the calculation of the variance of the tree lifetime we need the second 
moment of d which can be immediately obtained from the characteristic 
function. In the case of the model t we have 



ER 2 } 



duo 2 



6 2 



log 



y + i 
y (1 + 7). 



dy. 



J u=0 

By using some well known integral relations it can be shown that 

2 



ER 2 } = ~/3 2 7 



log 



i + i 

7 
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where Li 2 (u) = Y^k=\ u I s the so called Jonquiere's (dilogarithm) 
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Figure 8: Dependence of the variance and the relative variance of the tree lifetime -d on 
qx < 1 in models t and g. 

function. Taking into account this form of the second moment we can write 



D 2 W 



io K |i + i q2 



2 1-ga 



92 (1 - 9l) 



1^2 



1 92 

21-91 



(54) 
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Performing similar calculations in the case of the model g for the variance of 
■& g we obtain 



D ' 2 M = 1 " Ml - ^ - 14, (-j*-) . (55) 

In Fig. H one can see the dependence of the variance as well as the relative 
variance of the tree lifetime on the parameter qi in the cases of both models 
t and g. When qi is approaching to 1 from below the fluctuation of the tree 
lifetime becomes unlimitedly large, and so in the vicinity of the critical state 
the average lifetime loses almost completely its information content. 

4 Concluding remarks 

The probability distribution of the number of nodes with k > outgoing 
lines has been investigated in randomly evolving trees defined in Special 
attention was paid on the stochastic properties of end-nodes. We found that 
the birth and death of end-nodes in randomly evolving trees are playing a 
decisive role in the dynamics of the process. 

It is remarkable that the relative variance of the number of end-nodes vs. 
time has well-defined maximum when the evolution is either subcritical or 
supercritical. In the case of critical evolution the relative variance increases 
monotonously with the time. On the contrary, the relative dispersion of the 
number of dead end-nodes vs. time has a minimum just after the beginning 
of the evolution. The minimum can be seen in each of evolution states. 

We defined the lifetime of randomly evolving trees and derived a non- 
linear differential equation for the probability that the lifetime is larger than 
a given positive real number x. Three theorems have been proven to obtain 
asymptotic expressions for the survival probability of subcritical, critical and 
supercritical trees. Though the average lifetime of supercritical trees is always 
infinite it has been shown that the probability of finite lifetime of supercritical 
trees is larger than zero. In other words, randomly evolving supercritical trees 
may have finite lifetime. 
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